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We investigate properties of attractors for scalar field in the Lorentz violating scalar-vector-tensor 
theory of gravity. In this framework, both the effective coupling and potential functions determine 
the stabilities of the fixed points. In the model, we consider the constants of slope of the effective 
coupling and potential functions which lead to the quadratic effective coupling vector with the 
(inverse) power-law potential. For the case of purely scalar field, there are only two stable attractor 
solutions in the inflationary scenario. In the presence of a barotropic fluid, the fluid dominated 
solution is absent. We find two scaling solutions: the kinetic scaling solution and the scalar field 
scaling solutions. We show the stable attractors in regions of (7, £) parameter space and in phase 
plane plot for different qualitative evolutions. From the standard nucleosynthesis, we derive the 
constraints for the value of the coupling parameter. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

One of the interesting and widely exploited cosmologi- 
cal models is the scalar-tensor theory of gravity, the the- 
ory of a scalar field coupled to gravity. A remarkable 
phenomenon such as inflation, is usually discussed in the 
frame of this model. Recently, however, a scenario of 
implementing local Lorentz violation in a gravitational 
setting is considered to imagine the existence of a tensor 
field with a non- vanishing expectation value, and then to 
couple this tensor to gravity or matter fields. The vector 
field picks out a preferred frame at each point in space- 
time, and any matter fields coupled to it will experience 
a violation of local Lorentz invariance [H, 0. In a cos- 
mological background, such a vector field acts to rescale 
the effective value of Newton's constant ||. Moreover, 
from the study of the spontaneous breaking of Lorentz 
symmetry due to a vector field [3, PI H_0Kmany current 
experiments and observations @, M, USi Till . Qjf can be 
explained. In Ref. [l3[ it has been studied the late time 
evolution of the vector field perturbation and its effects 
on cosmological observables. They found that the CMB 
and matter power spectra are modified. 

More recently, a great interest has been devoted to 
study cosmological inflation in the framework of Lorentz 
violating scalar- vector-tensor theory of gravity [l4[ • They 
have shown that the Lorentz violating vector affects the 
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dynamics of the inflationary model. One of the interest- 
ing feature of this scenario, is the exact Lorentz violating 
inflationary solutions in the absence of the inflaton poten- 
tial. In this case, the inflation is completely associated 
with the Lorentz violation. Depending on the value of 
the coupling parameter, the three kinds of exact solu- 
tions are found: the power law inflation, de Sitter infla- 
tion, and the superinflation. Moreover, the dynamics of 
superinflation in the context of Loop Quantum Cosmol- 
ogy, in which the Friedmann equation is modified by the 
presence of inverse volume corrections has been studied, 
recently [l5[. 

The dynamical system of the scalar field with canoni- 
cal Lagrangian has been widely studied [l6|, [TtJ , among 
which the global structure of the phase plane has been 
investigated and various critical points and their phys- 
ical significances have been identified and manifested. 
General properties of attractors for scalar-field dark en- 
ergy scenarios which possess cosmological scaling solu- 
tions have been studied in Ref. [IH (see also 

Cosmological attractor solutions have been found and 
studied by several authors for various classes of poten- 
tials. The main purpose of this paper is to show that 
it is possible to find attractor solutions in the Lorentz 
violating scalar-vector-tensor models in which both the 
effective coupling function and the potential function are 
specified, and their relation matters. In other words, we 
will find a class of models in which the dynamics of the 
system dependents on the effective coupling function and 
the potential function, as well as on their relation. More 
specific, we consider the model in which the slope of the 
effective coupling vector and the potential in the Lorentz 
violation are constants. Then we find the Lorentz vio- 
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lating model with the quadratic effective coupling vector 
and the inverse or power-law-potential. The quintessence 
scenario with this potential has been well studied in the 
conventional universe [2(1 Hl[ . In the case of a tachyon 
field the corresponding potential is given by V = Vq4>~ 2 



j22Ll23. 24]. The dynamical attractor of the cosmological 
system has been employed to make the late time behav- 
iors of the model insensitive to the initial condition of 
the field and thus alleviates the fine tuning problem. In 
quintessence models, the dynamical system has tracking 
attractor that makes the quintessence evolves by tracking 
the equation of state of the background cosmological fluid 
so as to alleviating the fine tuning problem Qj], H(| • 
This paper is organized as follows. In Section[lll we set 
down the general formalism for the scalar-vector-tensor 
theory where the Lorentz symmetry is spontaneously bro- 
ken due to the unit-norm vector field. We derive the gov- 
erning equations of motion for the canonical Lagrangian 
of the scalar field. In Section IIII1 we study the attrac- 
tor solutions for the purely scalar field. In Section IIV1 
we extend our analysis in the presence of the barotropic 
fluid. The critical points of the system and their stability 
are presented. The final Section is devoted to the con- 
clusions. In the Appendix, we present the stability of the 
fixed point in which both the slope of effective coupling 
function and potential function are constant parameters. 



II. LORENTZ VIOLATING 
SCALAR- VECTOR-TENSOR 

In the present section, we develop the general re- 
construction scheme for the scalar-vector-tensor gravita- 
tional theory. We will consider the properties of gen- 
eral four-dimensional universe, i.e. the universe where 
the four-dimensional space-time is allowed to contain any 
non-gravitational degree of freedom in the framework of 
Lorentz violating scalar-tensor-vector theory of gravity. 
Let us assume that the Lorentz symmetry is sponta- 
neously broken by getting the expectation values of a 
vector field as < Olu^it^lO >= — 1. The action can be 
written as the sum of three distinct parts: 



S — Sq ~\~ S u ~\~ S(h 



(1) 



where the actions for the tensor field S g , the vector field 



S u , and the scalar field iS^, respectively, are given by 



16ttG 



(2) 



S,t, 



+A (u"tt p + 1)] , 
d x\/—g C<j, . 



(3) 
(4) 



In the above (i — 1,2,3,4) are arbitrary param- 

eters and C$ is the Lagrangian density for scalar field, 



expressed as a function of the metric g M „ and the scalar 
field 4>. A is a Lagrange multiplier. Then, the action (fl]) 
describes the scalar-vector-tensor theory. For the time- 
like vector field, we impose a constraint 



= -1. 



(5) 



Here, we take as the dimensionless vector, and ac- 
cordingly, (3i has the dimension of mass squared. Thus, 

1/2 

(3i 1 gives the mass scale of symmetry breakdown. The 
preferred frame determined by the vector differs from 
the CMB rest frame and the alignment of these frames 
had been achieved during the cosmic expansion as is ex- 
plained in the Appendix of Ref. [14]. In this setup, the 
preferred frame is selected through the constrained vector 
field and this leads to violating the Lorentz symmetry. 

For the background solutions, we use the homogeneity 
and isotropy of the universe spacetime 



ds 2 = -M 2 {t)dt 2 



e 2a( - t) 6 lj dx i dx j 



(0) 



where Af is a lapse function. The scale of the universe is 
determined by a. We take the constraint 



•" - ' ^,0,0,0 



(7) 



where Af = 1 is taken into account after the variation. 
Varying the action ((T|) with respect to g^, we have field 
equations 



(8) 



where T^ v = T^} + T$ is the total energy-momentum 

tensor, T$ and T$) are the energy-momentum tensors 
of vector and scalar fields, respectively, defined by the 
usual formulae 



(9) 



The time and space components of the total energy- 
momentum tensor are given by 



T u - -n 



Ti — Pu + V<l> 



(10) 



where the energy density and pressure of the vector field 
are given by 



Pu = -WH 2 , 

Pu 

f3 = /3 1 + 3/3 2 + ft 



3 + 2^ + 2^) iff' 



(11) 
(12) 

(13) 



Note that /J4 does not contribute to the background dy- 
namics. A prime denotes the derivative of any quantities 
X with respect to a. X' is then related to its derivative 
with respect to t by X' = (dXfdt)H- 1 = XH^ 1 where 
H = da/dt = a is the Hubble parameter. From Eqs. 
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(fTTj) and (fl2j) . one obtains the energy equation for the 
vector field u 



p'u + 3(p« - 
and for the scalar field 



Pv 



-3H 2 P' , 



(14) 



(15) 



The total energy equation in the presence of both the 
vector and the scalar fields is, accordingly, 



p' + 3(p + p) = 



P = Pv 



(16) 



Substituting Eq. ([10]) into the Einstein equations (JSJ), 
we obtain two independent equations, called the Fried- 
mann equations, as follows: 



1 



1 



8nG0 



8ttG/3 
(HH'-[ 



H 2 



H 2 ) = - 









i> ! 




' P<$> 






-H 2 










(17) 



3p<A 
J 



(18) 



The second term on RHS of Eq. (jlBf is a consequence 
of the coupling vector field as a function of scalar field. 
If 0i = 0, thus without the vector field, the above equa- 
tions reduce to the conventional ones. And in the case 
Pi = const., the above equations lead to the Friedmann 
equations given in Ref. [3j. 

Let us define the effective coupling function as follows 



= + 



1 



8nG 



then Eqs. (fTTj) and (fT5)) can be rewritten as 

1 



H z = — 

3/3 

W_ _ 



2 70 



(19) 

(20) 
(21) 



where the equation of state for the scalar field is defined 

by 70 = 1 +p<p/p^- 

For a given scalar field Lagrangian with the FRW back- 
ground, we can obtain the equations of motion for a scalar 
field by using Eq. (fTS"|) . Let us consider the Lagrangian 
density of a scalar field 4* with a potential V(4>) in Eq. JJJ : 



^ = -\{V4>? - V{ct>) , 



(22) 



where (V0) 2 = d^d v <\>. For the homogeneous field, 
the energy density p$ and the pressure of the scalar 
field may be found as follows 



V 

V 

1<I> 



60 



6/3 - (j) 1 ' 



= -2 



30 - <j>' 2 
60 - cj)' 2 



(23) 
(24) 
(25) 



Substituting Eq. (f23j) into Eq. (|2"0j) and also substituting 
Eq. (|25|) into Eq. (f2"Tjl . the Friedmann equations lead to 



(26) 
(27) 



H z 1 - = -^V 

6/3 / 3/3 



2/3 



Now differentiating Eq. (f2"3")) with respect to a and using 
Eq. (|15p . we obtain a dynamical equation for the scalar 
field 6, 



(28) 



which is subject to the Friedmann constraint given by 
Eq. 



ATTRACTOR SOLUTIONS FOR PURELY 
SCALAR FIELD 



III. 



Equations (|26| - (|28j) are the governing equations of mo- 
tion which we will use to study dynamical attractor for 
purely scalar field. We introduce the following dimen- 
sionless variables 1251: 



„2 _ 



Ai 



60 ' 

_ik 

s/0 

00 M 
~0\ 



A 2 



V 



m 2 



v 



VI 2 V )4> /V 



(29) 
(30) 

. (31) 



and, accordingly, the governing equations of motion could 
be reexpressed as the following system of equations: 



(32) 
(33) 



-3x(l-x 2 ) + x l-(X 1 + X 2 )y 2 , 



xy 



A', 



3X " V2 (Al + A ^ 
V6A 2 (v t -1\x 



X' 2 = -V6X2 ( r 2 - 1)» , 



(34) 
(35) 



where a prime denotes a derivative with respect to the 
logarithm of the scale factor, a = In a. 

In general, the parameters Ai, A2, Ti and T 2 are vari- 
ables dependent on <fr and completely associated with the 
Lorentz violation. In particular, Ai and Ti are purely 
Lorentz violation parameters and Ti can be written as 
a function of Ai in this case. By definition (|3"Tj) . T 2 can 
be written as a function of Ai and A2. Thus, in order 
to construct viable Lorentz violation model, we require 
that the coupling function and the potential function 
V should satisfy the condition Ti > 1/2 and T 2 > 1, 
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respectively. In this paper, we want to discuss the phase 
space, then we need certain constraints on the coupling 
function and potential function. In particular, we study 
the case of constants Ti and IV We take the form of Ti 
and T2 as rx = 1/2 and T2 = 1, respectively. Equations 
([34]) and (f35|) imply that Ai and A2 are nearly constants 
in this case. 

Integrating Eq. (f3"Tj) with respect to 0, we obtain a 
Lorentz violating model 

p(<j>) = ^ 2 , V(<f>) = Vo^ , (36) 

where £, n and Vo are parameters. Hence, 

Ai = -2^/1 , A 2 = ^Ai . (37) 



In the following, we study the case power-law potential 
n > 0. 

From the above equation, we require £ > 0. Then, the 
dynamical system of equations (|32[) ~ (|35[) can be rewritten 
as an autonomous system: 



x 1 = -3x(l-x 2 )-(n + 2)^m y 2 



xy . 



3x+(n + 2)\/^- 



The equation (|26[) leads to the constraint 
x 2 + y 2 = 1 . 

In term of the new variable, the equation (|27p reads 

— = -3x 2 - 2^6£x . 
rl 



(38) 
(39) 

(40) 
(41) 



Integrating this equation with respect to a one shows 
that all critical points, where x is a non-zero constant, 
correspond to an evolution of the Hubble parameter H 
given by H oc e _Q / p , where 



1 



P 



3x 2 + 2^Mx ' 



(42) 



An inflationary phase is one wherein the universe un- 
dergoes an accelerating expansion, i.e., the scale factor 
satisfies a > 0. Inflation ends when this condition is 
violated. In the frame of Lorentz violating scalar-vector- 
tensor, the condition for the accelerating universe is 



(43) 



where 7^ is an equation of state for the scalar field, 
-fcf, = 1+W0. When the condition (|4*3"|) is satisfied, we find 
three types of the inflationary solution given by Eq. (|42|) : 
the power-law inflation (p > 0), the de Sitter inflation 
(1/p = 0) and the superinflation (p — —\p\ < 0), de- 
pending on the values of £. We will show that there 



Power law 
inflation J 


de Sitter inflation 




/ Superinflation 



1 2 
Acceleration parameter 



FIG. 1: There are three possibilities of the inflationary solu- 
tions in the Lorentz violating scalar- vector-tensor theory: the 
power-law inflation, the de Sitter inflation and the superinfla- 
tion. 



exists a stable attractor for this three types of the infla- 
tionary solution in the framework of Lorentz violating 
scalar-vector-tensor theory of gravity. In the conven- 
tional case, £ — * 0, we are only left with a power-law 
inflation p — l/3a: 2 , where the condition for the accel- 
erated universe is, accordingly, 7^ < 2/3. In Fig. [TJ we 
show three possibilities of the inflationary solutions. Here 
we define an acceleration parameter, 



q = l 
that simplifies 



1 " ^ ~ 2V6£ 



P = 



l-q 



(44) 



(45) 



Then, the condition for the accelerating universe 
becomes q > 0. It can be shown that the inflationary 
(accelerated) solutions are classified into three cases (see 
Fig. [1]): a) < q < 1, a solution in which p > 0. In 
this case we have a power law inflation where a(t) ~ t p . 
b) q = 1, 1/p = 0. This is a de Sitter solution ait) ~ 
e Ht . c) q > 1. In this case, p = —\p\ < 0. Hence, 
the solution becomes a(t) ~ (— £)~' p ', t < 0. Thus, this 
solution represents the super-inflationary universe. 

In order to study the stability of the critical points, 
using the constraint Eq. (|66|) . we first reduce Eqs. (|38|) 
and (|39p to a one dimensional equation, 



x' = -3 



(I-* 2 ) 



(46) 



If we linearize the system about the critical points x — > 
x c + u, we could readily write the first-order perturbation 
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equation as 



2x\3x+(n + 2) 




3(1 -x 2 ) 



u , (47) 



which yields one eigenvalue. 

Depending on the values of £ and n, we find three 
critical points where x' and y' vanish: 

• Point (A±): 

{x Cl y c )=(±l, 0) is a kinetic dominated solution. 

• Point (B): 



(x c ,y c )=(-(n + 2)^/6, y/1 - (n + 2) 2 £/6) is a 
potential-kinetic dominated solution. 

Notice that the above critical points depend on the val- 
ues of the effective coupling parameter £ and the power 
of the potential n. For the positive values of £ and n, 
the linear perturbation shows that the point (A + ) is 
always unstable, whereas the point (A_) is stable for 
£ > 6/(n + 2) 2 . Thus, we only left with one possibility 
attractor solution in the kinetic dominated solution. The 
potential-kinetic dominated solution, point (B), is stable 
for < £ < 6/(n + 2) 2 . In this case, the point (A_) is 
unstable. 

For the point {A^), the universe accelerates for £ > 
1/6. Combining the stability of the critical point and the 
condition for the accelerating universe, the inflationary 
attractor solutions are given by < n < 4, £ > 6/(n + 2) 2 
and n > 4, £ > 1/6. Then, we find that the three types 
of inflationary (accelerated) solutions depend on the po- 
tential parameter and the effective coupling parameter. 



• The power-law inflationary attractor 
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2 < n < 4 
n > 4 , 



<«<! 



(n + 2) 2 
1/6 < £ < 3/8 . 



• The de Sitter inflationary attractor: 
n > 2 , 



• The super-inflationary attractor: 



< n < 2 , £ > 



6 



(n + 2) 2 ' 



n > 2 , 



(48) 
(49) 

(50) 

(51) 
(52) 



In the case potential-kinetic dominated solution, the 
accelerated expansions of the universe are given by < 
n < 2, £ > and n > 2, < £ < 2/(n 2 - 4). Then we 
find 



• The power-law inflationary attractor: 
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2<n<4 , 0<£< 



(53) 



(n + 2) 2 ' 

n>4 , 0<£<^-^. (54) 



• The de Sitter inflationary attractor: 

n = 2, 0<£<3/8. (55) 



• The super-inflationary attractor: 



0<n<2, 0<£< 



6 



(n + 2) 2 



(56) 



In summary, the effective coupling parameter and the 
potential model determine the stabilities of the critical 
points and the inflationary solutions. There are two pos- 
sibilities of attractor solutions given by the critical point 
(yl_) and (B). 



IV. ATTRACTOR SOLUTIONS IN THE 
PRESENCE OF BAROTROPIC FLUID 

For a realistic model, we consider the effect of an ad- 
ditional component. We carry out cosmological dynam- 
ics of the scalar field <fi in the presence of a barotropic 
fluid whose equation of state is given by p 7 = (7 — l)p 7 , 
where 7 is an adiabatic index, < 7 < 2. We assume 
that there is a barotropic fluid, not explicitly coupled 
to the scalar field and the vector field. Then, the to- 
tal energy-momentum tensor in Eq. ([5]) for this case is 



Ti m V f< 

LLU J- 111' 



(«) 



rr(<t>) 



■T$), where 



Tut* = {P~t +P~ l )n ll n v +p~ l g liv 



(57) 



is the energy-momentum tensor of the matter field. Here 
is the four velocity. The time and space components 
of the Einstein equations (|SJ) yield 



m 1 



8ttG (p u + p$ + p 7 ) , 



(58) 



and 



2H'H = -8nG (p u + p u + P4> + P4> + Pl + p 7 ) . (59) 



Substituting equations (TT|) and (|2"3"]) for the energy den- 
sity of vector field and scalar field into equation (|58[> . 
respectively, one finds 



H 2 = ± ( V + p 7 



3(3 \ 1 _ H 



(60) 



The evolution equation for a barotropic fluid is 

P'j = -37P7 ■ ( 61 ) 

The scalar field obeys the same equation of motion, 
Eq. (|28j) . The second Friedmann equation, Eq. ([59]) . be- 
comes 



EL 



2~J 



, Pi _ 
2H 2 fi 



(62) 



Equation (gS]) together with Eqs. ([gD ]) - ^ are the gov- 
erning equations of motion which will be used to study 
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dynamical attractor for a scalar field in the presence of 
the barotropic fluid. 

For the cases of the constants Ai and A2, the govern- 
ing equations can be written as the two-dimensional au- 
tonomous system: 



In the presence of barotropic fluid, the condition for 
the accelerating universe is 



^ x >\\j\ (7+^(7^-7)-- ) < 71) 



x! = -x [7(1 - x 2 - y 2 ) + 2x 2 ] - 3x 

+ ] jl[X 1 (l-x 2 )+X 2 y 2 ] , (63) 
y = \y [7(1 - x 2 - y 2 ) + 2x 2 ] 



^(X 1+ X 2 )xy. (64) 



The system of equations are symmetric under the reflec- 
tion [x, y) — + [x, —y). In what follows, we will restrict our 
discussion to the existence and stability of critical points 
to the upper half plane y > 0. In the case Pi =const., 
Ai — > 0, the scalar field dynamics in the Lorentz violat- 
ing scalar-vector-tensor theories is then reduced to the 
scalar field dynamics in the conventional one. But, the 
effective gravitational constant is rescaled by Eq. (fT9|) . 
In this case, the cosmological attractor solutions can be 
studied by a scalar exponential potential of the form 
V(<fi) — Vq exp(— \2<t>/-\fp) where j3 =const. This ex- 
ponential potential gives rise to scaling solutions for the 
scalar field [P?! ]. 

Also, the Friedmann constraint, Eq. (|60p . becomes 



P7 _ 
3H 2 f3 



(65) 



where the contribution of the scalar field to the total 
energy density is 



The equation (|T0[) can be simplified by 

1 



P 



l-q> 

where we have defined an acceleration parameter q. 



1 = 1 - Ifeff 



(72) 



(73) 



Thus the three types of inflation are also possible in the 
presence of the barotropic fluid. 



A. Stability of the fixed points 

The critical points (x c ,y c ) are obtained by imposing 
the conditions x' = and y' = 0. Substituting linear 
perturbation x — > x c + u and y — > y c + v about the critical 
points into Eqs. ([63| and ([64]) . we obtain, to first-order 
in the perturbation, the equations of motion 



U , = M U 
vl \ v 



(74) 



In what follows we clarify the properties of the five crit- 
ical points given in Table [T] (see Appendix) . We analyze 
the stability of the critical points with the background 
barotropic fluid < 7 < 2. The stable solutions of these 
critical points are 



a, 



m 2 i3 



2 1 2 

x +y 



In term of the new variable we find 
W 3 



H 



77e// 



(66) 



(67) 



where we have defined the effective equation of state of 
the universe: 



7e// = 7 + ^ (70 - 7) - 2y 3 Aix . (68) 

When x is a non zero constant, Eq. (|6"T)) corresponds to 
an evolution of the Hubble parameter given by 



H oc e~ a/p , 



(69) 



where 



a) Two kinetic dominated solutions 
For the point (A + ), the eigenvalues are 



mi = 3 + (n + 2)W y , m 2 = 3(2 - 7) + 2^ . (75) 

Thus, the critical point (A + ) is always unstable for n > 0. 
For the critical point (A-), the eigenvalues are 



mi = 3 - (n + 2)y y , m 2 = 3(2 - 7) - 2a/6£ . (76) 

The stable solutions for the kinetic dominated solution 
are 



< 7 < 



2n 
n + 2 ' 



3(2 - 7 ) 2 



(77) 



and 



p= 3 

27e// 



(70) 



2n 
n + 2 



<7<2 , 



(n + 2) 2 



(78) 
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In this kinetic dominated solution, we have the follow- 
ing relations: 



fy = l, 70 = 2, 7e// = 2 + 4 A /-^ , (79) 



(80) 



— = 1 - q A , q A = -2 + 2V6£ 

PA 



Note that the inflationary attractor solution dose not 
exist for 7 = 4/3. However, in the background matter 
(7 = I), the de Sitter inflationary attractor corresponds 
to the parameters, n > 3, < £. the superinflation- 
ary attractor does not exist in this case and the power- 
law inflationary attractor corresponds to the parameters, 
2 < n < 6, 3/4n < f < 3/8 and n > 6, 1/8 < £ < 3/8. 



b) Scalar field dominated solution, point (B) 
The eigenvalues are 

(n + 2) 2 

mi = -3+^ — ^£ > TO 2 = -3 7 + n(n + 2)£ . (81) 

The stable solutions for the scalar field dominated 
solution are 



0<7<^, 0<e<^^, (82) 



n + 2 



i(n + 2) ' 



and 



2™ „ . 6 
<7<2, < £ < — „, 2 



n + 2 



(n + 2) 2 



(83) 



The critical point exists for < £ < 6/(n + 2) 2 . In this 
case, we find the following relations: 

= 1, 7* = ^(n + 2) 2 £, 7e // = ^(n 2 - 4)£ (84) 



PB 



1 - 3B , 



l-^(n 2 -4)£ 



(85) 



The universe accelerates for 0<n<2,£>0 and n > 2, 
< £ < 2/(n 2 -4). 

c) Kinetic scaling solution, point (C) 
The eigenvalues are 

3 7 2< 3(2- 7 ) 4£ 

mi = y-2-^' m2 = — 2— + (2-^y (86) 

The stable solutions for the kinetic scaling solution are 
2n 



< 7 < 



n + 2 4n 
We obtain the following relations: 
8£ 



37(2-7) <e< 3(2-7)! . (87) 



3(2 - 7 ) 2 ' 



70 = 2 , Jeff = 7 - 



8£ 



3(2-7) ' 

3 7 , 4£ 



— = 1 - qc , gc = 1 - -7T , ,„ , 
PC 2 (2 - 7) 



(89) 
(90) 



The universe accelerates for 0<7<2/3,£>0 and 
2/3< 7 <2, £> (3 7 -2)(2- 7 )/8. 



d) Scalar field scaling solution, point (D) 
The eigenvalues are 



where 



T = — 3 1 — 



(n + 2)7 
2n 



(91) 



(92) 



D _ 3 [4< - 3 7 (2 - 7)] [3 7 - (n + 2)<] _ 
2n 2 £ 

The stable solutions for the scalar field scaling solution 
are given by 



< 7 < 



2n 
n + 2 



37 <c< 37( i: 7) . (94) 



n(n + 2) 



4n 



The critical point exists for < 7 < 2, < £ < 37(2 - 
7 )/4n. We find the following relations: 



(95) 



70 



PD 



37 2 



3 7 - 2< ' 



leff = , (96) 



1 - q D , Qd = 1 



n 

3 7 (n - 2) 
2n 



(97) 



and the universe accelerates for0<n<3,0<7<2 
and n > 3, < 7 < 2n/3(n - 2). 

From the above analyzes of the critical points, when 
one of the critical point is stable, then the other critical 
points are unstable or saddle depending on the values of 
7 and £. 



B. Phase-space diagrams 

In this subsection, we show phase-space diagrams of 
the stability solutions. Fig.[2]shows that different regions 
in the (7, £) parameter space lead to different qualitative 
evolutions for the potential model V(4>) = Vb^ 4 - 



Region A 



0<7<4/3, £>3(2- 7 ) : 



and 



4/3 < 7 < 2 



C>l/6 



(98) 



(99) 
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CI 
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1.0 
Y 



FIG. 2: Region of (7, £) parameter space for the potential 
model V(4>) = Vb0 4 . 



The kinetic dominated solution is stable in 

this case. The scalar field dominated solution and 
the kinetic scaling solution are saddle for £ > 1/6 
and £ > 3(2 — 7) 2 /8, respectively. The scalar field 
scaling solution is unstable. The effective equation 
of state is, accordingly, 7^ = 2 — 4(2£/3) x / 2 and 
the universe accelerates for £ > 1/6. For the in- 
flationary attractor solutions, we obtain: i) The 
power-law inflationary attractor solution (Region 
Al) corresponds to region of (7, £) parameter 
space: 1 < 7 < 4/3, 3(2 - 7 ) 2 /8 < ( , < 3/8 and 
4/3 < 7 < 2, 1/6 < £ < 3/8, ii) The de Sit- 
ter inflationary attractor solution (Line £ = 3/8) 
corresponds to region of (7, £) parameter space: 
l<7<2,£ = 3/8, and iii) The superinflation- 
ary attractor solution (Region A2) correspond to 
region of (7, £) parameter space: < 7 < 1, 
£ > 3(2 - 7) 2 /8 and 1 < 7 < 2, £ > 3/8. In 
Fig. [3J we show the phase plane plot for 7 = 1 
and £ = 1. The superinflationary attractor is the 
kinetic dominated solution. 



• Region B 



< 7 < 4/3 , < £ < ^ , 



and 



4/3< 7 <2, 0<£<l/6 



(100) 



(101) 



The scalar field dominated solution is stable in this 
case and it exists for £ < 1/6. The kinetic domi- 
nated solution and the scalar field scaling solution 
are saddle for < £ < 1/6 and < £ < 7/8, re- 
spectively. The kinetic scaling solution is unstable. 
The universe accelerates for < £ < 1/6. We only 
find the power-law inflationary attractor solution 



in this case, and it corresponds to region of (7, £) 
parameter space: < 7 < 4/3, < £ < 7/8 and 
4/3 < 7 < 2, < £ < 1/6. In Fig. gj we show 
the phase plane plot for 7 = 1 and £ = 1/16. The 
power-law inflationary attractor is the scalar field 
dominated solution. 



• Region C 



< 7 < 



3 7 (2- 



16 



^<£<^- 



7) 2 



(102) 



The kinetic scaling solution is stable in this case. 
The kinetic dominated solution, the scalar field 
dominated solution, and the scalar field scaling so- 
lution are saddle for £ < 3(2 - 7 ) 2 /8, £ > 1/6, 
and 37(2 — 7)/ 16 < £, respectively. Also, for the 
inflationary solutions, there are three cases to be 
considered, i.e., i) The power-law inflationary at- 
tractor solution (Region CI) corresponds to re- 
gion of (7, £) parameter space: < 7 < 1, 
37(2 - 7 )/16 < £ < 37(2 - 7 )/8 and 1 < 7 < 4/3, 
37(2-7)/16 < £ < 3(2- 7 ) 2 /8, ii) The de Sitter in- 
flationary attractor solution (Line £ = 37(2 — 7)/8) 
corresponds to region of (7, £) parameter space: 
< 7 < 1, £ = 37(2 -7) /8, and iii) The superinfla- 
tionary attractor solution (Region C2) corresponds 
to region of (7, £) parameter space: < 7 < 1, 
3 7 (2 - 7 )/8 < £ < 3(2 - 7 ) 2 /8. Fig. [5] shows the 
phase plane plot for 7=1 and £ = 1/4. The power- 
law inflationary attractor is the kinetic scaling so- 
lution. 

• Region D 

The scalar field scaling solution exists for < 7 < 
2, < £ < 37(2 - 7)/16, and it is stable for 

o< 7 <±, m 

The kinetic dominated solution is unstable. The 
scalar field dominated solution and the kinetic scal- 
ing solution are saddle for £ > 7/8 and £ < 
37(2 — 7)/6, respectively. The universe accelerates 
for < 7 < 4/3. For the inflationary attractor 
solutions, we only find the power-law inflationary 
attractor solution in this case, corresponding to the 
region of (7, £) parameter space: < 7 < 4/3, 
7/8 < £ < 37(2 - 7)/16. In Fig. EJ we show the 
phase plane plot for 7 = 1 and £ = 2.5/16. The 
power-law inflationary attractor is the scalar field 
scaling solution. 



C. Cosmological implications 

As we have shown in the previous section, there are 
some possibilities that lead to an accelerated expansion: 
the kinetic dominated solutions and the scalar field dom- 
inated solutions together with scaling solutions. The 



FIG. 3: The phase plane for the potential model V = Vb</> 4 , 
7 = 1, and £ = 1. The attractor is a kinetic dominated 
solution with x = —1, y = 0, where Qj, = 1, 7^ = 2, 7,;// = 
2(1 - 2^/2/3) and g = 2(^6 - 1). 



FIG. 5: The phase plane for the potential model V = Vo4> 4 , 
7 = 1, and £ = 1/4. The attractor is a kinetic scaling solution 
with x = — \/2/3, y = 0, where = 2/3, 7^ = 2, 7 e // = 1/3 
and g = 0.5. 





FIG. 4: The phase plane for the potential model V = Vb</> 4 , 
7 = 1, and £ = 1/16. The attractor is a scalar field dominated 
solution with x = — \/3/8, y = \/5/8, where fi^ = 1, 7^ = 
3/4, 7e// = 1/4 and g = 5/8. 



FIG. 6: The phase plane for the potential model V = Vb</> 4 , 
7=1, and £ = 2.5/16. The scalar field dominated solution 
is the saddle at x = — \/l5/4, y = 1/4, and the attractor is 
a scalar field scaling solution with x = — \/3/5, y = l/\/l0, 
where fi^ = 7/10, 7^ = 12/7, 7 e // = 1/2 and g = 1/4. 



later, however, is not viable to explain an accelerated 
universe at the present epoch since ~ 0.65 ± 0.05 
with ujfj, < —0.4 [27]]. In this section we consider the cos- 
mological implications of the scaling solutions. For the 
case power- law potential, we find two scaling solutions: 
the kinetic scaling solution and the scalar field scaling 
solution. 

In the background radiation, we find the following re- 
lations for the kinetic scaling solution: 



^0 = 6£ , 70 = 2 , j eff 



so 



30 . (104) 



The universe accelerates for £ > 1/6. However, this scal- 
ing solution is not viable to explain the accelerated uni- 
verse in the background radiation, because the universe 
is unstable in this case. 

The scalar field scaling solution has the density param- 
eter, 



"■■-siel'-T 



(105) 



for 7 = 4/3, where ranges £ between < £ < 2/3n. We 
also obtain 



70 



3(2 - <) 



Jeff 



4(n - 2) 
3n 



(106) 



Moreover, using the standard nucleosynthesis and the 
observed abundances of primordial nuclides, the strong 
constraint is that the fractional energy density of scalar 
field at temperatures near 1 MeV is < fl™ ax . This 
implies 



n(nn™ ax + 2) 
For SVV; ax ~ 0.2, the nucleosynthesis bound requires 
20 



(107) 



n{n + 10) < ^ ' 
for the Lorentz violation to be relevant. 



(108) 
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-10 12 3 

In a 

FIG. 7: The evolutions of the effective equation of state as a 
function of In a for the case £ = 1/32, 7 = 1 and n = 4. We 
choose the initial conditions = 0.01, y% = 0.12. 

From Eq. (fl~5)) , we define 7^ ^ = 74, + H 2 (3' / which 
is related with the effective equation of state parameter 
ui^^ by the relationship = + H 2 f3' / p^. In 

Fig. [7] we show the evolutions of equation of state as 
a function of a = In a for the case £ = 1/32, 7 = 1 
and n = 4. We choose the initial conditions Xi = 0.01, 
2/j = 0.12. The effective equation of state for the scalar 
field can cross —1. 



(7, £) parameter space. Based on this parameter space, 
we have shown the phase plane plot for the four differ- 
ent qualitative evolutions, i.e. the kinetic dominated so- 
lution, the scalar field (potential-kinetic) dominated so- 
lution, the kinetic scaling solution and the scalar field 
scaling solution. Another result is that the scalar-vector- 
tensor theory of gravity is still relevant at the nucleosyn- 
thesis. From this data, we derive the constraints for the 
value of the coupling parameter in Eqs. (|107|) . 

More interestingly, the model (|36[) allows for the 
inverse-power potential, n < 0. As an example, in the 
case n = —2, we only have one attractor solution in the 
potential dominated solution. The effective equation of 
state parameter 7^ is of the cosmological constant. Two 
attractor solutions, the scalar field dominated solution 
together with the scaling solution exist for n < — 2. 

When Ai is dynamically changing quantity, the effec- 
tive coupling vector could track the evolution of the scalar 
field. There exists one stable critical point that gives an 
acceleration of the universe at late time with the equa- 
tion of state parameter u> < —1/3. Thus, we can have a 
dark energy scenario in the frame of scalar-vector-tensor 
theory of gravity (28|. Finally, we would like to empha- 
size that there exists an attractor solution in the Lorentz 
violating scalar-vector-tensor theory of gravity. 



V. CONCLUSIONS 



In this paper, we have studied the attractor solutions of 
the scalar field in the frame of Lorentz violating scalar- 
vector-tensor theory of gravity. In this model, because 
of the dynamics of the effective coupling vector, the dy- 
namics of the scalar field is modified from the conven- 
tional cosmology. The characteristic of the slope, p, and 
the condition of the accelerating universe are given by 
Eqs. (HO]) and l[7IJ). respectively. 

For the model (|36[) without a barotropic fluid, the crit- 
ical points of a dynamical system completely depend on 
the coupling parameter with typical potential. There ex- 
ist two kinetic dominated solutions and one potential- 
kinetic dominated solution. Different regions in the 
parameter space lead to different qualitative evolutions. 
There are three typical inflationary solutions which are 
obtained from Eq. (|4"2")l: the power-law inflation, the de 
Sitter inflation and the superinflation. Three types of in- 
flationary attractor exist both in the kinetic dominated 
solution and scalar field dominated solution. Note that 
the supcrinflationary is stable in the case purely Lorentz 
violation without potential, where £ > 3/2. 

We have extended the dynamical system of the evolu- 
tion to a realistic universe model with a barotropic fluid. 
The stability of the critical point was discussed in Sec. 
IV. In this case, the fluid dominated solution is absent, 
instead, we find two scaling solutions: the kinetic scaling 
solution and the scalar field scaling solution. We have 
completely have shown the stable attractor in regions of 
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APPENDIX: STABILITY OF THE CRITICAL 
POINTS FOR CONSTANTS Ai AND A 2 

Let us consider the cases Y\ = 1/2 and T2 = 1. From 
Eqs. (f3"4"|) -([3l) )) . we obtain Ai and A2 nearly constants. 
The critical points (x c ,y c ) are obtained by imposing the 
conditions x' = and y' = 0. Substituting linear pertur- 
bation x —f x c +u and y — > y c +v about the critical points 
into Eqs. (Rl3"|) and ([M]) . to first-order in the perturbation, 
one obtains the equations of motion 

00 -*(:)• f*- 1 ' 
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TABLE I: The properties of critical points with barotropic fluid in the Lorentz violating scalar-vector-tensor theory of gravity 
for constants Ai and A2. 



Point 


X 


y 


Existence 


stability 


(A+) 


+ 1 





All Ai, A2 and 7 


Stable: Ai > 6 - VGX'z/VG and Ai > 3(2 - j)/V& 


(A-) 


-1 





All Ai, A2 and 7 


Stable: Ai < -(6 + V^X 2 )/V^ and Ai < -3(2 - *y)/VG 


(B) 


^(Ai + A 2 ) 




" (A!+A 2 ) 2 " 
1 6 


1/2 


(Ai + A 2 ) 2 < 6 


Stable: See Eq. (|A.12| 


(C) 







All A 2 , 7 2 and Ai / 


Stable: See Eci. (|A.15|I 


(D) 


(3/2) 1/2 ^ 




"37(2-7) 

2A 2 . A 2 


1/2 


A1A2 < 


Stable: See Eas. (1A.17|I-(|A.18|) 



TABLE II: The scalar field density parameter the equation of state 7^, the effective equation of state 7e//, the slope p and 
the acceleration parameter, q, for the cases of constants Ai and A2. 



Point 


Q4, 


l<t> 


Jeff 


1/p 


q 


(A+) 


1 


2 


2 - 2^2/3Xi 


3- v^Ai 


-2 + y^Ai 


(A-) 


1 


2 


2 + 2^2/3Ai 


3 + v^Ai 


-2 - v^Ai 


(B) 


1 


i (Ai + A2) 2 


HAS -A?) 


i(AS-A?) 


1-KA^-Af) 


(C) 


2\'{ 


2 


2A^ 




1 3 7+ K 


3(2-7) a 


7 3(2-7) 


2 ' (2-7) 


2™ T (2-7) 


(D) 


37 ( I AxA 2 \ 


7(l-^) _1 






i-M 1 -^) 



where the elements of the matrix M are 
Mn 



M12 
Mai 
M 22 



-3 + 3(2 - 7 )x 2 + - [7 (1 - x 2 

— v^Aix , 
(\/6A 2 - 372:^ y , 



3(2- 7 )x-\/-(A 1+ A 2 ) 



2/ ■ 



-372/ 2 + g [7 (1 ~ ^ - y 2 ) + 2^ 2 



(Ai +A 2 )x 



2a; 2 J 
(A.2) 
(A.3) 

(A.4) 



(A.5) 



The critical points together with the stability analy- 
sis for constants Ai and A 2 are shown in Table HI For 
completeness, the scalar field density parameter f2^, the 
equation of state 7^, the effective equation of state 7 e //, 
the slope p, and the acceleration parameter q are shown 
in Table [II] The eigenvalues of the stability matrix M, 
Eq. (jA.ip , are as follows: 

• Point (A±) 

(x c ,y c ) = (±1,0) is a kinetic dominated solution, 



(A.6) 
(A.7) 

(A.8) 



mi = 3 =F y 2 ( A i + ^2) 
m 2 = 3(2 - 7) T \/6Ai . 
The fixed point (A + ) is stable for 



Ai > 



6- V6A 2 



and Ai > 



3(2-7) 
V6 



and the fixed point (A-) is stable for 



6±\/6A 2 3(2-7) 
Ai < 7= , and Ai < ,y . (A.9) 



x/6 



V6 



• Point (B) 

(x c ,y c )=((\i + A 2 )/V6, [1 - (Ai + A 2 ) 2 /6] 1/2 ) is a 
scalar field dominated solution, 



1 2 

mi = -3+2 ( Al + A 2) ' 



(A.10) 

m 2 = -3 7 + A 2 (Ai +A 2 ) . (A.ll) 



The fixed point is a stable for 



— \/6 — A 2 < Ai < \/6 - A 2 , and 
Ai > 37 ^ 2 for A 2 < , or 



Ai < 



37 -A 2 
A 2 



for A 2 > 



(A.12) 



• Point (C) 

(x c , y c )=((2/3) 1 / 2 Ai/(2 — 7), 0) is akinetic scaling 
solution, 



3 /„ \ An 

mi = --(2- 7 ) + 



m 2 



2 V " 2-7 
37 AiA 2 

T ~ 2-7 ' 



(A.13) 
(A.14) 
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The fixed point is stable for 
3 



The fixed point is a stable for 



A?>^(2- 7 ) 2 , 



and 



Ai > 7 -* for A 2 < , or 



Ai < 



2A 2 

37(2 - 7) 
2A 2 



for A 2 > 



(A.15) 



2-7 
Ai < J -\ 2 , 

7 

37(2-7) 
2A, 



Ai < 



• Point (D) 

(x c ,y c )=((3/2)VVA 2 ,[37(2-7)/2Ai-A 1 /A 2 ] 1/2 ) where A 2 < 0, 



is a scalar field scaling solution, 
3 



mi 2 



1± 



\ 



i 8(3 7 2 - 67 + 2AiA 2 )(3 7 -AiA 2 -Ai) 



3A1 



2-7(1+^ 



2-7 
Ai < ^-A 2 , 

7 

37(2-7) 



Ai > 



2A 2 



(A. 16) where A 2 > 0. 



and 



Ai > 



37 z3 

A 2 



and 

Ai < 



37 -A| 



(A.17) 



(A.18) 
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